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ABSTRACT 

It is possible to visualize the Cosmic Web as an interconnected network of one- 
dimensional filaments, two-dimensional sheets and three-dimensional volume-filling 
structures which we refer to as clusters.We have used the Local Dimension D, which 
takes values D = 1,2 and 3 for filaments, sheets and clusters, respectively, to anal- 
yse the Cosmic Web in a three-dimensional volume-hmited galaxy sample from the 
Sloan Digital Sky Survey Data Release 7. The analysis was carried out separately 
using three different ranges of length-scales: 0.5-5, 1-10 and 5-50 h^^Mpc. We find 
that there is a progressive increase in the D values as we move to larger length- 
scales. At the smallest length-scale, the galaxies predominantly reside in filaments 
and sheets. There is a shift from filaments to sheets and clusters at larger scales. Fila- 
ments are completely absent at the largest length-scale (550 h^^Mpc). Considering 
the effect of the density environment on the Cosmic Web, we find that the filaments 
preferentially inhabit regions with a lower density environment as compared to sheets 
and clusters which prefer relatively higher density environments. A similar length- 
scale dependence and environment dependence was also found in a galaxy sample 
drawn from the Millennium Simulation which was analysed in exactly the same way 
as the actual data. 
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1 INTRODUCTION 

Galaxy redshift surveys provide us with a picture of the large scale structures in the present day 
universe. All the major galaxy redshift surveys like the Center for Astrophysics (CfA) Survey 
(Geller & Huchra, 1989), the Las Campanas Redshift Survey (LCRS) (Shectman et al, 1996), the 
Two-Degree Field Galaxy Redshift Survey (2dFGRS) (CoUess et al. 2001) and the Sloan Digital 
Sky Survey (SDSS) (York et al. 2000) clearly show that the galaxies are distributed in a complex 
interconnected network of filaments, sheets and clusters encircling nearly empty voids. This inter- 
connected network is often referred to as the "Cosmic Web" (Bond et al., 1996). Quantifying the 
Cosmic Web and understanding it's origin is one of the most interesting and challenging issues in 
cosmology. 

A wide variety of statistical measures have so far been employed to quantify the Cosmic Web. 
The void probability function (White, 1979), the percolation analysis (Shandarin & Zeldovich, 
1983) and the genus curve (Gott, Mellot & Dickinson, 1986) are some of the earliest statistics 
introduced to quantify the topology of the galaxy distribution. The Minkowski functionals (Mecke 
et al., 1994) provide a global characterization of structures. Ratios of the Minkowski function- 
als can be used to define the 'Shapefinders' which are a set of shape diagnostics for both simple 
and topologically complex objects (Sahni et al. 1998). Bharadwaj et al. (2000) have introduced 
a two dimensional (2D) version of Shapefinders. There are quite a few different techniques that 
have been introduced to identify voids in the galaxy distribution (El- Ad & Piran, 1997; Hoyle & 
Vogeley, 2002; Platen et al., 2007; Neyrinck, 2008). Stoica et al. (2007) have proposed a three 
dimensional object point process to delineate filaments in the large scale structures. Colombi, 
Pogosyan & Souradeep (2000) have proposed a Hessian based statistics to study the topology of 
excursion sets at the percolation threshold. The Smoothed Hessian Major Axis Filament Finder 
(SHMAFF) (Bond et al. 2010) has been introduced to identify filaments, sheets and clusters in the 
galaxy distribution. Aragon-Calvo et al. (2007) have introduced a multiscale morphology filter to 
automatically segment the cosmic structures into its basic components namely clusters, filaments 
and walls. Sousbie et al. (2008) have proposed a skeleton formalism to quantify filamentary struc- 
tures in a three dimensional density field. Aragon-Calvo et al. (2008) have introduced the Spine of 
the Cosmic Web which provides a complete framework for the identification of different morpho- 
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logical components of the Cosmic Web. Each of the techniques mentioned above quantifies one 
or almost a few aspects of the complex network referred to as the Cosmic Web. A comprehensive 
quantification of the Cosmic Web is still forthcoming, leaving considerable scope for work in this 
direction. 

The early galaxy redshift surveys like the LCRS probed thin, nearly two dimensional (2D) 
slices through the universe. Initial investigations (Bharadwaj et al., 2000) which used the 2D 
Shapefinders show the galaxy distribution in the LCRS to have excess filamentarity in comparison 
to a random distribution of points. Bharadwaj et al. (2004) show the filamentarity in the LCRS 
to be statistically significant up to a length-scales of 70 — 80 Mpc but not beyond. Longer 
filaments, though present in the data, were not found to be statistically significant, and were pos- 
sibly the outcome of chance alignments of shorter filaments. Bharadwaj & Pandey (2004) have 
used N-body simulations to show that the filamentarity in the LCRS is consistent with the LCDM 
model with a mild bias. Subsequent work (Pandey & Bharadwaj, 2005) performed a 2D analysis 
of thin slices from the SDSS. This confirmed the results obtained earlier using the LCRS. It was 
also found that the distribution of brighter galaxies has a lower connectivity and filamentarity as 
compared to the fainter ones. The filamentarity was also found to depend on other galaxy proper- 
ties like colour and morphology (Pandey & Bharadwaj, 2006), and the star formation rate (Pandey 
& Bharadwaj, 2008). Pandey & Bharadwaj (2007) have computed the filamentarity for galaxy 
samples in different luminosity bins and compared these to the filamentarity in simulated galaxy 
samples with different values of the linear bias parameter to obtain a luminosity-bias relation. A 
recent 2D analysis of the luminous red galaxies in the SDSS shows the filamentarity to be statis- 
tically significant to length-scales as large as 100 to 130 h"^ Mpc (Pandey et al., 2011), which are 
considerably larger than those found earlier in the LCRS and the SDSS main galaxy sample. 

Though the Cosmic Web is relatively easy to visualize and analyze in 2D, this has several 
limitations which can only be overcome in a three dimensional (3D) analysis. For example, a 2D 
filament could actually be a section through a 3D sheet. Sarkar & Bharadwaj (2009)(henceforth 
Paper I) have proposed the Local Dimension as a 3D statistics for analysing the Cosmic Web. 
This can, in principle, be used to classify different structural elements along the Cosmic Web as 
filaments, sheets and clusters. Tests with cosmological N-body simulations (Paper I) show that 
the structures identified a filaments by the Local Dimension match quite well with the visual 
appearance. The structures identified as sheets, however, could not be visually identified. This was 
attributed to the fact that the visual appearance is determined by the most dominant structures in 
the field which usually are the filament, whereas the sheets which are relatively diffuse structures 
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are not visually identified. Paper I also showed that the Local Dimension could also be used to 
address a variety of issues like determining the relative fraction of galaxies in filaments, sheets 
and clusters respectively. 

In the present work we have used the Local Dimension to analyse the patterns in the galaxy dis- 
tribution in the Sloan Digital Sky Survey Data Release Seven (SDSS DRV). We have also carried 
out a similar analysis on a galaxy catalogue from the semi-analytic model of galaxy formation im- 
plemented in the Millennium Simulation (Springel et al., 2005), and compared the results from the 
actual data with those from the simulation. The Local Dimension classifies the different structural 
elements along the Cosmic Web as filaments, sheets and clusters. We have used this to study the 
relative fraction of galaxies in these three different kinds of structures. Further, we also study how 
these three different kinds of structures are distributed with reference to varying density environ- 
ments. For example, it is possible that the filaments are preferentially distributed in high density 
environments relative to the low density environments. It is well accepted that the Cosmic Web is 
a complex network whose morphology and connectivity will depend on the length-scale at which 
the analysis is carried out. The galaxy distribution is expected to approach homogeneity at large 
length scales, and it is known (Sarkar et al., 2009) that the SDSS DR6 is consistent with homo- 
geneity at length-scales beyond 60 — 70 h"^ Mpc. Given these considerations, we have carried out 
the entire analysis for three different ranges of length-scales namely 0.5 to 5 /i~^Mpc, 1 to 10 
h-^Mpc and 5 to 50 h-^Mpc. 

An brief outline of the paper follows. In Section 2 we describe the Local Dimension which is 
the method that we have adopted to analyse the Cosmic Web and in Section 3 we discussed the 
data and the method of analysis while Section 4 contains the results and conclusions. 

2 THE LOCAL DIMENSION 

The Cosmic Web may be thought of as an interconnected network of different structural elements. 
Any particular structural element may be classified as being either a cluster, a filament or a sheet. 
The Local Dimension, introduced in Paper I, is a simple yet effective method to quantify the 
shape of individual structural elements along the Cosmic Web. The entire analysis is based on 
the following argument. Consider, for example, a galaxy "G" located in a filament which is , by 
definition, a one dimensional structure. We expect A^(< R), the number of other galaxies within a 
sphere of comoving radius R centered on 'G', to scale as A^(< R) = AR. Similarly, we expect a 
scaling A'^(< R) — AR^ with D — 2 and 3 if G were located in a sheet and a cluster respectively. 



Exploring the Cosmic Web in SDSS DR7 5 
The exponent D quantifies the dimension of the galaxy distribution in the neighbourhood of G, 
and hence we refer to it as the Local Dimension. In principle the Local Dimension provides a 
technique to determine whether a galaxy is located in a filament, a sheet or a cluster. 

Consider, for example, a part of the Cosmic Web where there is a filament connected to a sheet 
(Figure 1 of Paper 1). We expect D — 1 and 2 for the galaxies located in the filament and the sheet 
respectively. We also expect that the galaxies located near the junction of the two structures will 
not exhibit a well defined scaling behaviour, and it will therefore not be possible to determine the 
Local Dimension D for these galaxies. In a typical situation where the Cosmic Web is a complex, 
interconnected network of different kinds of structural elements only a fraction of all the galaxies 
in the entire survey will have a definite value of the Local Dimension D. 

We have, until now, not taken into account the fact that galaxies are discrete objects with finite, 
non-zero intergalactic separations. It is thus necessary to consider a finite range of comoving length 
scales i? < i?2 in order to determine if there is a scaling behaviour A^(< R) = AR^ and thereby 
estimate the Local Dimension D. Contrast this with a continuum where it is possible to study the 
scaling behaviour in the limit R ^ and it is not necessary to refer to a finite range of length- 
scales. Further, the Cosmic Web is expected to look different when viewed at different length- 
scales. We also expect the Cosmic Web to approach a homogeneous network beyond length-scales 
of 60 — 70 Mpc (Sarkar et al., 2009) where the galaxy distribution is known to approach 
homogeneity. It is thus interesting and useful to separately analyze the scaling behaviour across 
different ranges of length-scale. Based on this we also introduce a lower length-scale Ri and study 
the scaling behaviour over different ranges of length- scales Ri < R < R2. 

The values of Local Dimensions once determined, by fitting power law over the length-scales 
Ri < R < R2, can be put to use to address a variety of issues. First, we can identify individual 
structures like filaments or sheets. The relative abundance of different D values allows us to esti- 
mate the fraction of galaxies that reside in sheets, filaments and clusters respectively. The spatial 
distribution of the D values allowed us to study how the different kinds of structural elements are 
interconnected or "woven into the Cosmic Web". 

3 DATA AND METHOD OF ANALYSIS 
3.1 SDSS DR7 Data 

Our present analysis is based on galaxy redshift data from the SDSS DRV (Abazajian et al., 2009). 
The SDSS DR7 includes 11,663 square degrees of imaging and 9380 square degrees of spec- 
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troscopy with 930, 000 galaxy redshifts. For the present work we have used the Main Galaxy Sam- 
ple for which the target selection algorithm is detailed in Strauss et al. (2002). The Main Galaxy 
Sample comprises of galaxies brighter than a limiting r band Petrosian magnitude 17.77. The data 
was downloaded from the Catalog Archive Server (CAS) of SDSS DRV using a Structured Query 
Language (SQL) search. We have identified a contiguous region in the Northern Galactic Cap 
which spans —40° < A < 33° and —30° < rj < 30° where A and rj are survey co-ordinates 
defined in Stoughton et al. (2002). A volume limited galaxy subsample was constructed in this 
region by restricting the extinction corrected Petrosian r band apparent magnitude to the range 
14.5 < rrir < 17.77 and restricting the absolute magnitude to the range —19 < Mr < —20.5. 
This gives us 47705 galaxies in the redshift range 0.035 < z < 0.076 which corresponds to the 
comoving radial distance range 104 < r < 223 Mpc. At the redshift z — 0.035, the co- 
moving volume corresponding to our sub-sample subtends 124 /i~^Mpc x 104 /i~^Mpc along the 
transverse directions. 

3.2 Millennium Data 

The Millennium Simulation (Springel et al., 2005) is a large cosmological N-body simulation 
which traced 2160^ particles from redshift 127 to the present in a periodic comoving box of side 
500 h~^Mpc. We have used a semi-analytic galaxy catalogue generated by Guo et al. (2011) who 
updated the previously available galaxy formation models (Springel et al. 2005; Croton et al. 2006; 
De Lucia & Blaizot 2007) with improved versions and implemented galaxy formation models on 
the Millennium Simulation. Semi-analytic models are simplified simulations of the formation and 
evolution of galaxies in a hierarchical clustering scenario incorporating all relevant physics of 
galaxy formation processes. The spectra and magnitude of the model galaxies were computed 
using population synthesis models of Bruzual & Chariot (2003). Using the peculiar velocities, we 
map the galaxies to redshift space and then identify a region having the same geometry as our 
actual data. Applying the same magnitude cuts as those used for the actual data, we have extracted 
the same number of galaxies as in our final SDSS data and used this in our subsequent analysis. 

3.3 Method of analysis 

We have determined A^(< R) at several R values for each galaxy in our sample. For each galaxy, 
its distance from the survey boundary sets the largest value of R for which it is possible to estimate 
N{< R). We have assigned the Poisson error AN{< R) — ^yN{< R) to each measured value 
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Figure 1. The galaxy distribution witliin 5 h^^ Mpc of one of the SDSS DR7 galaxies that was identified to have D 1 across the length-scales 
0.5 to 5 /i~^Mpc. The galaxies are all shown with crosses. We have used a Friend-of-Friend (FoF) algorithm with a linking length of \/^h~^ Mpc 
to highhght any connected structure in the galaxy distribution in this figure. This yields a single filamentary structure containing the majority of the 
galaxies. All the galaxies belonging to this connected structure are shown with cells. The filled cells represents galaxies with D « 1. It was not 
possible to determine D for the other galaxies in this figure. 



of A^(< R). A minimisation procedure was used to determine the best fit power law A^(< 
R) = AR^ to the A'^(< R) measured for each galaxy. The value of D is accepted as the Local 
Dimension corresponding to the particular galaxy if the chi-square per degree of freedom of the 
power law fit satisfies, / < 1- The power law fit is rejected for larger values of / and the 
Local Dimension is undetermined for these galaxies. The fitting procedure was restricted to values 
within the range i?i < < i?2 in length-scale. 

We have carried out the analysis for three different ranges of length-scales, each covering a 
decade. Values of D were determined separately across the length-scales 0.5 to 5 ^~^Mpc, 1 to 10 
/i~^Mpc and 5 to 50 /i^^Mpc. It was possible to determine a definite value of D for 3484, 9082 
and 288 galaxies at the three respective ranges of length-scale mentioned above. 

In order to illustrate our method of classifying cosmological structures we show the galaxy 
distribution in the vicinity of one of the SDSS galaxies that was identified to have D ^ 1 across 
the length-scales 0.5 to 5 /i~^Mpc (Figure 1). The expected filament is clearly visible passing 
through the center of the figure. We note that there are several galaxies with D = 1 arranged along 
the filament visible in the figure. Figure 2 shows a similar plot for one of the galaxies with D Ri 1 
from the Millennium data. 
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4 RESULTS AND CONCLUSIONS 
4.1 Distribution of D values 

We first analyse the fraction of centers witli different D values, shown in Figure 3. The D val- 
ues were divided into bins of width ±0.25. The error bars in the data have been estimated using 
bootstrap re-sampling of the data. Ten bootstrap samples were used for this purpose. 

The solid curve in Figure 3, which corresponds to the results for 0.5 to 5 h~^Mpc, shows a 
broad peak with a maxima at D — 1.5. The bin centered at D = 2 contains the second largest 
fraction of galaxies. These two bins together contain more than 50% of the centers for which D 
could be determined. This value indicates that the galaxies in the Cosmic Web are predominantly 
contained in sheets and filaments at the length-scales 0.5 to 5 /i^^Mpc, with the sheets being 
somewhat more dominant than the filaments. The dashed curve in Figure 3, which corresponds 



Exploring the Cosmic Web in SDSS DR7 9 




to length-scales 1 to 10 /i~^Mpc, shows a sharp peak at D = 2. The fraction of galaxies in the 
neighbouring bins (D = 1.5 and 2.5) falls to nearly half the fraction in this bin. The three bins 
at D = 1.5, 2 and 2.5 together contains more than 80% of the centers for which D could be 
determined. This indicates that the galaxies in the Cosmic Web are predominantly in sheets over 
the length-scale 1 to 10 h~^Mpc. The dot-dashed curve in Figure 3, which corresponds to the 
length-scales 5 to 50 h~^Mpc, peaks at D — 2.5 and D — 3. The D — 2.5 contains the maximum 
no of centers. The two bins combinely contains more than 70% of the centers. This indicates that 
the galaxies in the Cosmic Web are predominantly in sheets and clusters (volume filling structures) 
over the range of length- scales 5 to 50 h~^Mpc. Further, it is interesting to note that we do not find 
any center with D — lov 1.5 at this range of length-scales. This indicates the complete absence of 
filamentary structures at the largest range of length-scales (5 to 50 h'^Mpc) that we have probed. 

There is a shift to larger D values in Figure 3 as we progressively consider larger length- 
scales. The entire curve showing the fraction of centers as a function of D shifts to the right 
when we consider larger length-scales. It is quite evident from this that the nature of the structural 
elements that make up the Cosmic Web differs depending on the length-scale at which we view 
the Cosmic Web. At small scales (0.5 — 5 /i~^Mpc) we have a mixture of sheets and filaments. The 
fraction of sheets increases at larger length-scales (1 — 10/i~^Mpc), whereas we predominantly 
have a combination of sheets and clusters at the largest length-scale (5 — 50 h~^Mpc) that we have 
probed. Filaments are completely absent at the largest length-scale. 

For comparison, we have also applied the Local Dimension to analyse the galaxy distribution 
in a semi analytic galaxy catalogue from the Millennium Simulation. We have used the semi 
analytic galaxy catalogue (Quo et al., 2011) to extract three different data samples with exactly 



1 Sarkar P. , Pandey B. , Bharadwaj, S. 



'■ SDSS 


' 1 ' 1 




1 ' 1 ' ■ 

D=1 

- — D=2 
D=3 




























1 ' 


1 . 1 







2 4 6 8 10 



Figure 5. The three curves which correspond to D = 1 , 2 and 3 respectively show the fraction of galaxies as a function of u. The results are for 
the SDSS data using Ri = 0.5 h'^ Mpc, R2 = 5h-^ Mpc and = 1.58 h'^ Mpc. 
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Figure 6. The three curves which correspond to D = 1, 2 and 3 respectively show the fraction of galaxies as a function of v. The results are for 
the SDSS data using Ri = 1 h'^ Mpc, R2 = 10h-^ Mpc and Rs = 3.16 fe-^ Mpc. 



the same geometry and galaxy number density as our SDSS sample. These three simulated data 
samples were analysed in exactly the same way as the actual data. The results showing the fraction 
of centers at different D values are presented in Figure 4. The three different curves in this figure 
correspond to the same range of length-scales as three different curves in Figure 3. We find that 
the fraction of centers with different D values have very similar distributions in the actual SDSS 
data and the Millennium Simulation. These curves for the SDSS data are nearly identical to those 
from the Millennium Simulation except for the fact that at largest length-scale (5 — 50 h^^ Mpc), 
the fraction of centers with D = 2.5 is larger in the SDSS data as compared with the Millennium 
simulation. This difference in the fraction of centers lies in the 1 — a error-bar. 
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Figure 7. The three curves which correspond to D = 1 , 2 and 3 respectively show the fraction of galaxies as a function of u. The results are for 
the SDSS data using Ri=5h-^ Mpc, R2 = 50 /i"! Mpc and Rs = 15 /i"! Mpc. 



4.2 Environment dependence 

In the previous subsection, we liave seen that the nature of the structural elements that make up the 
Cosmic Web changes depending on the length-scale at which we view the web. Our investigations 
show that there is a progressive transition from filaments to sheets and then clusters as we go 
from smaller to larger scales. In this subsection, for a fixed range of length-scales, we investigate 
if the distribution of a particular kind of structural element is related to the density environment. 
As mentioned earlier, the nature of the structural elements (ie. D value) depends on the range of 
length-scale used to determine D. The density environment too depends on the length-scale at 
which we smooth the density field. The length-scale range Ri to R2 used for estimating D should 
thus be consistent with Eg, the smoothing length-scale for the density field. In our analysis we have 
chosen Eg to be the geometric mean of Ei and E2 i.e Eg — \fE[E^. The geometric mean was 
chosen instead of the algebraic mean because the latter is expected to be more biased towards E2. 
We expect the geometric mean to give a more representative estimate of the range of length-scales 
El to E2 which span a decade. 

For fixed values of Ei, E2 and Eg we focus on the distribution of the centers with a particular 
D value. For this purpose, the D values were divided into bins of width ±0.5 centered dXD — 1,2 
and 3. To determine the density environment, we have converted the entire galaxy distribution to a 
density field on a grid of spacing [0.5 /i~^Mpc]^ using the Cloud-in-Cell method. This density field 
is then smoothed with a Gaussian kernel having a smoothing length Es- The smoothing is carried 
out in Fourier space by multiplying the Fourier transform of the density field with exp(— A;^i?^/2) 
and transforming back to real space. The density field at any grid point was quantified using the 
dimensionless ratio v — 5/a where 5 — is the density contrast of the smoothed density field 
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at the particular grid point and a is the standard deviation of the density contrast of the smoothed 
density field evaluated using all the grid points that lie within the survey volume. Considering only 
the centers for which it is possible to determine a D value, we use the value of u at the grid point 
nearest to the center to assign a u value to each of these centers. The value of u associated with 
any of the galaxies gives an estimate of the density environment in the vicinity of the structural 
element centered on that galaxy. 

We first consider the range, Ri = 0.5 Mpc and _R2 = 5 Mpc for which i?^ = 1.58/i^^Mpc 
and a = 5.15. Figure 5 shows the results for this range of length-scales. We first consider only 
the centers with D = 1 for which the results are shown in the solid curve of this figure. This 
curve shows the fraction of centers with a particular u value. The u values were divided into bins 
spanning ±0.25 for this purpose. The dashed and dash-dotted curves respectively show the cor- 
responding results for D = 2 and 3 respectively. We find that the distribution of the fraction of 
centers as a function of v is qualitatively similar for all three values of D. The fraction shows a 
peak near ~ 1, with a very sharp decline in the fraction at i/ < 1 and a relatively gradual decline 
at u > 1. While the behaviour is qualitatively similar, there are quantitative differences between 
the different D values. We see that the fraction peaks at a somewhat smaller u value fovD — l 
as compared to D = 2 and 3. Atu < 1 the values of the fraction are somewhat larger for D = 1 
as compared to D = 2 and D = 3. Further, the values of the fraction in the vicinity of the peak 
are somewhat larger for D = 2 as compared to D = 3. The behaviour is reversed for u > 1. The 
fraction is smaller for = 1 in comparison to D = 2 and 3 in the range 1 < u < 3. For u > 3, 
we find that D — 3is roughly consistent with D — 1, whereas the fraction is considerably higher 
for D = 2. 

Figure 6 shows the same quantities as Figure 5 with the difference that the range of length- 
scales now corresponds to i?i = 1 and i?2 = 10 h~^Mpc for which Rg = 3.16 Mpc and 
a = 3.09. The behaviour, we find, is very similar to that in Figure 5 except that the peak has 
shifted to a smaller u value (u ~ 0.5). For z/ < 0.5, the fraction is larger for D = 1 in comparison 
to D — 2 and 3. This is reversed for i/ > 0.5 where the fraction is smaller fovD — 1 in comparison 
to D = 2 and 3. There is another transition around i/ ~ 2, where for D = 3 the fraction falls below 
that of D — 2. We find that, for i/ > 2, D = 1 and 3 are comparable and are lower than D — 2. 

Figure 7 shows the result for i?i = 5 and i?2 = 50 Mpc for which Rg = lb.d>h~^ Mpc and 
a = 0.97. The u values were divided into bins spanning ±0.5 for this purpose. The behaviour is 
different compared to that in Figure 5 and Figure 6. The analyse for this length-scale shows total 
absence of center in a filament. For v < 1.5, the fraction for D = 3 larger in comparison to D — 2. 
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Figure 8. The three curves which correspond to D = 1 , 2 and 3 respectively show the fraction of galaxies as a function of v. The results are for 
the data from the Millennium Simulation using _Ri = 0.5 h'^ Mpc, R2 = 5 /i"! Mpc and Rs = 1.58 h'^ Mpc. 
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This situation is reverse for u > 1.5, where the fraction is smaller for D = 2 in comparison to 
D = 3. 

We may interpret the different curves in Figures 5 and 6 as representing the probability of 
finding a particular kind of structural element in the density environment corresponding to u. For 
example, the curve for D=l in Figure 5 gives the probability of finding a filament in the interval 
1/ — 0.25 to z/ + 0.25 for different values of u. Similarly, the curves for D=2 and D=3 show the prob- 
ability of finding a sheet and a cluster respectively. Our analysis shows that the filaments, sheets 
and clusters have different probability distributions. The filaments are preferentially distributed in 
low density environments relative to the distribution of sheets and clusters. We have a cross-over 
from this behaviour at intermediate densities. The sheets are preferentially distributed relative to 
filaments and clusters in the high density environments. The u value where these transitions occur 
depends on our choice of Ri and i?2- For 0.5 to 5 h"^ Mpc, the density contrasts ranges are u < 1, 
1 < u < 3 and u > 3 while for 1 to 10 Mpc the u ranges are u < 0.5, 0.5 < u < 2 and 
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Figure 10. The three curves which correspond to D = 1 , 2 and 3 respectively show the fraction of galaxies as a function of f. The results are for 
the data from the Millennium Simulation using = 5/1-1 Mpc, ^2 = 50/1-1 Mpc and i?s ^ 15.8/1-1 Mpc. 



u > 2. These findings indicate that the way in which different structural elements are distributed 
along the Cosmic Web depends jointly on two factors (a) the local density environment, and (b) 
the length-scale at which we view the Cosmic Web. 

For comparison, the above analysis was also performed using the galaxy distribution in the 
Millennium Simulation. The analysis was carried out in exactly the same way as for the actual 
SDSS DRV data. The results for the Millennium simulation shown in Figures 8, 9 and 10 are 
analogous to the plots shown for the SDSS DR7 in Figures 5, 6 and 7 respectively. We find that 
the results from the Millennium Simulation are very similar to those obtained for the actual SDSS 
DR7 for length-scale 0.5 to 5 h^^Mpc and 1 to 10 /i~^Mpc. The results for the length-scale 5 to 
50 h~^Mpc from Millennium Simulation differs to those obtained from the SDSS DR7. The v 
values, in Figure 10, were divided into bins spanning ±0.5. At u < —0.5 the values of the fraction 
are somewhat similar for D = 2 and D — 3. The fraction is smaller for D = 2 as compared to 
D = Sin the range —0.5 < u < 1. This trend reversed for i/ > 1, where D = 2 is larger compared 
to D = 3. 

In summary the analysis of both the SDSS DR7 and the Millennium Simulation exhibit similar 
trends. The Local Dimensions were estimated separately in three different ranges of length-scales, 
0.5 to 5 h-^Mpc, 1 to 10 h-^Mpc and 5 to 50 /i~^Mpc. We find that there is a progressive shift 
in the D values as we move to larger length-scales. At small length-scales there is a mixtures of 
sheets and filaments, and the fraction of sheets increases as we move to larger length-scales (1 to 
10 /i~^Mpc). We find that sheets and clusters are the predominant structures at the largest length 
-scales (5 to 50 /i~^Mpc). Filaments are absent at this length-scales. It is interesting to note that 
Forero-Romero et al. (2009) find that the mass fraction of filaments decreases while that of sheet 
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increases with an increase in the smoothing length-scale, which is consistent with the results of 
this paper. 

The gradual transition, with increasing length-scales, from filaments {D — 1) to sheets (D — 
2) and then to clusters (D — 3) is indicative of and consistent with a transition to homogeneity 
(D — 3) at large length-scales. A recent analysis of the SDSS Main galaxy sample (Sarkar et 
al., 2009) shows that there is transition to homogeneity at 60 — 70 Mpc, and that the galaxy 
distribution is consistent with a homogeneous point distribution at length-scales larger than this. 
An earlier study (Pandey & Bharadwaj, 2005) had analysed thin, nearly 2-D, sections through 
the galaxy distribution and had found evidence for statistically significant filamentary patterns 
to length-scales as large as 80 Mpc. This is apparently inconsistent with the findings of the 
present paper which fails to find any filaments that span across the length-scale 5 to 50 Mpc. It 
should however be noted that the structures which were identified as filaments in the 2-D sections 
are quite likely to be sheets when the structures are viewed in full 3-D. 
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